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. Abstract. The status of the theoretical effort to calculate charge-symmetry-breaking pion produc- 

tion reactions will be presented. The main emphasis will be on the dd — > an reaction, especially 
in relation to the recent IUCF experiment. Suggestions for future theoretical and experimental work 
will be made. 



INTRODUCTION 

> 

This presentation is based mainly on [1], work done in collaboration with C. J. Horowitz, 
A. Nogga, A. C. Fonseca, C. Hanhart, G. A. Miller, J. A. Niskanen, and U. van Kolck. 

The Lagrangian of quantum chromodynamics (QCD) is almost symmetric under the 
exchange of the up and down quarks, which is called charge symmetry (CS). However, 
since the quark masses are different, this symmetry is broken, charge symmetry breaking 
(CSB). For hadrons, CS implies the invariance of the strong interaction under, e.g., the 
exchange of the proton and neutron. 

Experimental evidence for CSB has been demonstrated (before 2002) in p°-co mixing, 
the neutron-proton and other hadron mass splittings, the binding-energy difference of 
mirror nuclei such as 3 H and 3 He, the different scattering lengths of elastic nn and pp 
scattering, and in the difference between the proton and neutron analyzing powers of 
elastic np scattering. A review about CSB and related issues can be found in Ref._[2]. 

The n-p mass difference can be split into two parts, m n — m p = 8M + 8M 
1 .29 MeV/c 2 , where 8M is due to the strong interaction and stems from the up-down 
quark mass difference, while 8M is due to electromagnetic effects. The latter is known 
to be negative and various models gives typically dM ~ — 1 MeV/c 2 [3], but its precise 
value is not known. How can we determine these constants? The Lagrangian of chiral 
perturbation theory (#PT) contains two terms parametrized by 8M and 8M: 

^csb = -^nUt 3 -^-^)N+—NUt 3 + J N. (1) 



O 



o 

=5 



X 



VI ) 2 V VI 

The pion terms may give a different combination of 8M and 8M than for m n — m p , 
which makes it possible to use reactions involving pions to extract 8M and 8M. A first 
suitable choice is the CSB forward-backward asymmetry of np — > dn°, where originally 
a negative value Aft, = —0.28% was predicted from r\-iz Q mixing [4]. When the 8M 
and 8M terms were included the sign changed and Aft, = 0.23%-0.60% [3]. This latter 




FIGURE 1. Feynman diagram for leading order pion production in the dd — > an reaction (simplified 
model), indicating the labeling of nucleons and defining basic kinematic variables. The cross indicates the 
occurrence of CSB and the dot a leading-order CS vertex. 



prediction is of the same sign and magnitude as the subsequent experimental result from 
TRIUMF: A ft = [17.2±8(stat) ±5.5(sys)] x 1(T 4 [5]. 

Another possiblity is the dd — > an reaction, which is CSB since the deuteron and a 
particle are self-conjugate under CS, while the pion wave function changes sign. Thus 
this reaction could not occur if CS were conserved. The cross section is proportional to 
the square of the CS amplitude, which is unique since all other CSB observables involve 
interferences with CS amplitudes. A long series of experimental work, culminating with 
the efforts at the Saturne accelerator at Saclay, lead only to decreasingly small upper lim- 
its [6] . Finally, the Saturne group claimed to see dd — > an at T d = 1 . 1 GeV [7] . This was 
refuted by members of the same collaboration who argued that the dd — > ayy reaction 
was responsible for the signal [8]. The importance of this background was later con- 
firmed by calculation of the double-radiative capture [9], using a model based on a very 
successful treatment of dd — > ann at similar energies [10]. Thus, the claimed signal is 
most likely a misinterpretation of a heavily-cut smooth dd — > ayy background [9]. This 
realization influenced the design of the IUCF experiment, which last year reported a very 
convincing dd — ► 0C7t° signal near threshold (o = 12.7 ± 2.2 pb at T d = 228.5 MeV and 
15.1 ±3.1 pb at 231.8 MeV), superimposed on a smooth dd — > ayy background [11]. 
This background is roughly a factor two larger than calculations based on Ref. [9], but 
has the expected shape. The data are consistent with the pion being produced in an s- 
wave, as expected from the proximity of the threshold (T d = 225.6 MeV). 

It is hoped that these two experimental results, when analyzed in a consistent theo- 
retical framework, can be used to extract the values of 8M and 8M [3]. Our aim here is 
to provide the first study of CSB in the near threshold dd — > an reaction using chiral 
techniques. Previous models of dd — > Ct7t° have been at higher energies and have used 
T]-7t mixing only [12, 13]. We are guided in our quest by general symmetry considera- 
tions, relying on chiral power counting to give a list of possible amplitudes, and trusting 
that realistic wave functions can be constructed for the four-body states. 



POWER COUNTING 



The relevant scales are the typical momentum for pion production {% = p/M ~ 
\Jm % lM) ordering the strong diagrams and the fine structure constant (a 1/137) 
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FIGURE 2. Left: NNLO diagrams with strong CSB [(a)-(k)]. Encircled vertices are sub-leading. Not 
all possible time orderings are shown. Right: Some typical higher-order diagrams with strong CSB [(1) 
and (m)]. A double line represents a A-isobar. Diagram (1) appears at N 3 LO whereas diagram (m) is a 
N 4 LO contribution. 



ordering electromagnetic effects. Starting with the strong diagrams, at leading order 
(LO) there is only the % rescattering diagram of Fig. 1, which is driven by 8M and 8M. 

At next-to-next-to-leading order (NNLO) there are several diagrams as displayed in 
Fig. 2. The first (a) (one-body) term is parametrized by 7]-K° mixing (resonance satura- 
tion) jSi = g ri {K Q \H\r\) /m^, unfortunately with large uncertainties in and (n°\H\ri). 
The pion exchanges (b)-(d) are smaller. We have not yet included the pion loops (e)-(k). 

The N*LO diagram of Fig. 2 are in #PT given by low-energy constants (LECs) that 
should eventually be fitted to data. Since the LECs are not known we chose to model 
diagram (b) by assuming resonance saturation, i.e., heavy-meson (a, p, co) exchanges. 
Also p-co mixing and hard photon exchanges contribute. This is very similar to the 
situation in pp — > ppK° where short-range mechanisms were also needed [14, 15]. The 
CSB parameter is again jSi . The A-excitations (1) has not been calculated. 

The electromagnetic leading order causes CSB in the wave functions. Higher orders 
contain CSB pion production. These amplitudes have not yet been included in our model. 

SIMPLIFIED MODEL 

As a first step toward a full model we chose to use simplified (Gaussian) wave functions, 
without initial state distortions, to gain insight into the problem [1]. The diagram for the 
simplified model is given in Fig. 1. The deuterons need to be slowed down (the cm. 
momentum is p ~ 460 MeV/c) in order to produce a pion and a particle at rest. Thus 
we would like to have momentum sharing between the deuterons so that it does not push 



the nucleons far out in the high-momentum tail of the a -particle wave function. 

With spatially symmetric wave functions the spin-isospin parts are anti-symmetrized: 



\d n ) = (12)i[12] , |J 34 ) = (34) 1 [34]o 
l«) = ^{((12) 1 (34) 1 )o[[12]o[34]o]o-((12)o(34)o)o[[12] 1 [34] 1 ] } 

where (i,j) s ([i,j]r) are the spin (isospin) Clebsch-Gordan couplings, with magnetic 
quantum numbers suppressed, for nucleons, or nucleon pairs, i and j coupling to spin 
s (isospin T). Amplitudes easily match \dd) with the first part of \a), while the second 
term is matched only if both spin and isospin of both nucleon pairs are flipped simultane- 
ously. Thus symmetries decide about matching, but do they support momentum sharing? 

The symmetries forbid LO pion exchange between the deuterons and only the NNLO 
recoil piece (—m 7t /M)p N survives [Fig. 2(c)]. Thus the LO is strongly suppressed. 

The NNLO one-body cannot allow for momentum sharing and only its recoil part 
can contribute to s-wave production, but it is supported by the symmetries and adds 
coherently over all nucleons. It is quite small, of the order 1 pb. 

The short range N 4 LO amplitudes survive the symmetries, allow for momentum 
sharing, and add coherently over all nucleons and with each other, dominating the cross 
section in our simplified model, as demonstrated in Table 1 . 

RESULTS 

The diagrams calculated so far are summarized in Table 1, where a comparison is also 
made with the plane wave impulse approximation (PWIA) using realistic wave func- 
tions, (CD-Bonn+Urbana or AV18+Urbana). The pion rescattering recovers somewhat, 
since new symmetries possible with realistic wave functions allow for momentum shar- 
ing. The one-body term increases dramatically since it does not have any momentum 
sharing and the high-p tail of the a wave function is no longer a tiny Gaussian. 

While the simplified model gets quite close to the experimental value, the realistic 
wave functions, without distortion, overshoot considerably. There are several reasons 
why this may happen. Firstly, this can be because the initial state interactions are not 
included yet. Secondly, the ignored diagrams (e.g., the pion loops at NNLO) may cause 
destructive interferences. Thirdly, since the strength is mainly driven by the one-body 
term, the resonance saturation procedure may not be accurate? To test this we turned 
off the r\-K mixing, which gave more reasonable cross sections. Another concern is the 
difference between the CD-Bonn and AVI 8 calculations. However, our dependence on 
potential choices can be assessed only when the full model is in place. 

CONCLUSION AND OUTLOOK 

We find that, at least in our simplified model, the LO diagrams are suppressed, while 
subleading amplitudes are important and dominate the cross section. Thus, these short- 
range terms may be important for np — > dn° and should be evaluated there as well. 



TABLE 1. Cross sections evaluated for dd — > GOT at the lowest IUCF 
energy, for simplified wave functions (SWF) and realistic (PWIA), 
either CD-Bonn or AVI 8. The experimental cross section is also given. 



Operator 
(mech.) 


SWF 
[pb] 


PWIA (CD-Bonn) 
[pb] 


PWIA (AV18) 
[pb] 


7i(8,8) 


0.011 


1.619 


1.384 


l(ri-n) 


0.688 


12.561 


10.341 


a(ri-7z) 


0.187 


0.826 


0.467 


C0(T]-7t) 


0.404 


1.084 


0.759 


p<J)-K) 


0.081 


0.185 


0.092 


p-co 


1.645 


4.191 


2.647 


y(EM) 


1.486 






CSC l(EMin wf) 




0.095 


0.080 


total 


22.955 


79.832 


46.508 


no 77-7T 


1.925 


9.065 


3.049 


Expt. 


12.7 ±2.2 



We are working on the completion of our model, especially regarding the dd initial 
state. Polarized dd — > dd data currently being analyzed at IUCF will be used to test the 
scattering wave functions. The remaining diagrams, the NNLO pion loops, photon two- 
and three-body terms, and possible A (p-wave pion) contributions, will be included. The 
dependence on the choice of wave function and potentials can then be tested. 

With realistic bound state wave functions the cross section is overestimated, which 
may reflect our ignorance of the unknown LECs. To address this issue, new experiment 
are necessary, e.g., np — > dft and unpolarized or polarized dd — > an , with angular 
distributions and possibly at higher energies. As a dd — > an experiment is being 
planned at COSY, hopefully new exciting data will be available within a few years. 

This work was supported in part by the DOE grants DE-FG02-93ER40756 and DE- 
FG02-02ER41218. 
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